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We study the emission of gravitons by a homogeneous brane with the Gauss-Bonnet term into 
an Anti de Sitter five dimensional bulk spacetime. It is found that the graviton emission depends 
on the curvature scale and the Gauss-Bonnnet coupling and that the amount of emission generally 
decreases. Therefore nucleosynthesis constraints are easier to satisfy by including the Gauss-Bonnet 
term. 
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I. INTRODUCTION 

In recent years there has been considerable interest 
in the suggestion that our universe is a brane: a sub- 
space embedded in a higher-dimensional bulk spacetime. 
In these models, ordinary matter is confined to our 
brane while the gravitational field propagates through 
the whole spacetime. Of particular importance is the 
Randall-Sundrum (RS) model, where a single brane is 
embedded in an infinitely extended AdS^ spacetime [if. 
At low energies, the zero mode of the 5D graviton is lo- 
calized on the brane, because of the strong curvature of 
the bulk due to a negative bulk cosmological constant, 
and 4D gravity is recovered. 

A natural extension of the RS model is to include higher 
order curvature invariants in the bulk. Such terms arise 
in the Ads/CFT correspondence as next-to- leading or- 
der corrections to the CFT Particularly, in the 
heterotic string effective action, the Gauss-Bonnet (GB) 
term arises as the leading order quantum corrections to 
gravity. This gives the most general action with second- 
order field equations in five dimensions [^, 0| and is in- 
vestigated in areas such as black holed and brane- world. 
The graviton is localized in the GB brane-world Q and 
deviations from Newton's law at low energies are less pro- 
nounced than in the RS model ■ 

Brane cosmologies with and without GB term has been 
investigated [aM E3] • Due to the cosmological symme- 
tries, most GB brane-world scenarios assume that the 5D 
spacetime metric is the generalized Schwartzschild-Anti 
de Sitter (Sch-Ads), described by the metric [12|. Il3| 
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ds 2 = -g(r)dt 2 + —— + r 2 fl lj dx' l dx : > , 
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where fiy is the three dimensional metric of space with 
constant curvature k = — 1,0,1, A is the bulk cosmolog- 
ical constant, a > represents the GB coupling. In the 
a — > limit, this reduces to the usual Sch-Ads metric. If 



the bulk is empty then C is necessarily constant in time. 
However particle interactions can produce gravitons that 
are emitted into the bulk at high energies on a brane. 
Therefore, in a realistic cosmological scenario, there ex- 
ists an avoidable bulk component and so C is no longer 
constant. This problem has recently been studied for a 
RS brane 0, [l^ . In this paper we examine the radiating 
GB brane-world and find what effects including the GB 
term has on the evolution of C. 

The rest of this paper is organized as follows: in section II 
we derive the energy loss through graviton radiation; in 
section III wc derive the emission rate of the bulk gravi- 
tons; in section IV wc numerically solve the system of 
equations under some approximations; in section V some 
conclusions are drawn. 



II. THE BULK AND THE BRANE 

In order to model the bulk spacetime metric, we use 
the five dimensional and Gauss-Bonnet generalization of 
the Vaidya metric given by [f3 |. 
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where v = const are ingoing plane-formed null rays. If C 
does not depend on v, then the metric (2) is a rewriting 
of the generalized Sch-Ads metric (1), as can be seen by 
the coordinate transformation v = t + J dr/f{r). From 
now on we assume that the brane is outside the horizon 
(/ > 0) and that the brane universe is spatially flat. 
The Vaidya type metric is a solution to Einstein-Gauss- 
Bonnet equations 



where 



Gab + 2aH a b + Ag a b = K T a b, 



Gab — Rab — ^ffob-R, 
H a b = RRab — %RacR C b ~ 2R cd R a 



(3) 



(4) 



i d cde d 
+ K a Kbcde 



-^gabLGB, 



(5) 



and the bulk energy-momentum tensor has null-radiation 
form, 



T ci b = ipk a kb- 



(6) 



Here, k 2 = 1/M 3 is the five dimensional gravitational 
coupling, ip is, for a brane observer, the flux of gravitons 
leaving a radiation dominated brane and k a is a null vec- 
tor. Thus, in our model the bulk gravitons are presumed 
to be emitted radially. By inserting the metric (2) and 
the stress energy tensor (6) into Einstein-Gauss-Bonnct 
equations (3) we find the evolution equation for C: 



dC 
dv 



2K 2 V>r 3 



h h 



The appropriate normalization of k a is given by k a u a = 
— 1, where u a is the brane's velocity vector. This implies 
that the only nonvanishing component is k r = k v = 1/v, 
where v — dv/dr and r is cosmic proper time on the 
brane. From u a u„ = 1 we obtain 



fv = r + s/f 2 + f. 



(8) 



In order to determine the behavior of the brane we 
have to impose the gen eralized Israel junction conditions 
which are given by [l5(, 

[A^] - + 2a(3LV] - h^[J] - 2P wva \KP a \) 

= -K 2 S^,($) 



where 



J ab = X -(2KK ac K\ + K cd K cd K a 



-2K nr K cd K 
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Here, K ab is the extrinsic curvature, hp V is the induced 
metric on the brane and Sp U is the brane energy momen- 
tum tensor. From these junction conditions we obtain 
the following Friedmann equation 12j, |15[ : 



c + + c--2 
8a 



(12) 



where 
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and cr represents the energy density of the matter source. 
The requirement that the standard form of the Fried- 
mann equation is recovered at sufficiently low energy 
scales leads to the relation 



2 
^4 



1 



M 2 Pl (1 + 7)*' 



(14) 



where Mpi is the reduced 4D Planck scale, i~ 2 = (1 — 
sjl + 4aA/3)/4a is the AdS curvature scale, 7 = 4a/* 2 , 
and we have the standard assumption that the energy 
density on the brane is separated two parts, the ordinary 
matter component, p, and the brane tension, A, such that 
(j = p + A. We also assume zero cosmological constant on 
the brane, 



k 2 \ 



2(3-7) 



(15) 



The GB term is considered as the lowest-order stringy 
correction to the 5D Einstein action, so the GB energy 
scale should be larger than the RS energy scale. From 
this consideration, we have 7 < 0.15 [l6| . 
Then the Raychaudhuri equation is written as 



H + H = - 
2 ( 6 1 / 3 - 



1 - b 1 / 3 2C n 2 j> 
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(16) 



where b 1 / 3 = y/l + 4aA/3 + 8aC/r 4 . 

Because of graviton emission, the brane energy is not 

conserved, 



p + 4:Hp 



(17) 



The factor of 2 on the right hand side is due to the fact 
that the brane is radiating a flux of gravitons into both 
sides. 



III. PRODUCTION RATE OF BULK 
GRAVITONS 

In order to determine quantitatively the energy loss ip 
we follow the same procedure as in the RS case [9(. First, 
we evaluate the cross section of the process <f> + <^> — ► KK 
graviton, where is a particle confined on the brane. To 
compute this cross section we have to check whether or 
not the cosmological influence is negligible. In the GB 
high energy regime the Hubble rate is approximated as 



2 



H ~ T 4 / 3 /« 1/3 M, where T is the temperature of the 
brane particle. Here, the GB energy scale should be 
smaller than M so that there is the GB regime before 
the quantum gravity regime. From this condition we 
have a ^S> M~ 2 16j. Therefore, the temperature T is 



bigger than the Hubble rate H in the GB regime if we 
assume T -C M. Even after the GB regime T ^> H as 
shown in Thus, we find that the cosmological influ- 
ence can always be neglected. 

Let us consider the linear perturbations of the GB metric, 

ds 2 = e - 2A(z) {(rj^ + h^dafdx" + dz 2 } , (18) 

in axial gauge with 4d TT condition. Here, A(z) = 
log(|z|/£ +1). We decompose the graviton into KK 
modes, 



hfiv = j dmu m (z)4>nv(x), (19) 
where the modes u m (z) are given by jl7j |. 



Urn (z) 



A 



lm(\z\+i) 



2(1 



[Y 2 (m(\z\+£)) 



A 2 ) 

+AJ 2 (m(\z\ +£))] 
Y x (m£) + xmlY 2 (m£) 
' JJmt) +xm^J 2 (m^)' 



(20) 
(21) 



and satisfy the normalization J dzu^ n (z)u m i (z) = 8(m — 
ml). Y and J are the Bessel functions and Neumann 
functions respectively, and x = 7/(1 — In the a — > 
limit we recover the RS result. The coupling of the bulk 
graviton to the brane matter is described by the action 



S int = k(1 - 7 ) 1/2 / dmu m (0) / d A xS^4> 



(22) 



The overall factor (1 - 7) 1 / 2 is a GB correction @, E3- 
From this action we can calculate the amplitude for the 
scattering of brane particles leading to a KK emission. 
This calculation is quite analogous to the procedure al- 
ready described in the context of flat extra dimensions 
[ill, the only difference being the coupling constant in 
(22). Using those results, one finds that the spin and 
particle-anti particle averaged squared amplitude is given 
by 



Y>\M\' 



K 2 a 



7 )\u m (0)\ 2 A-, 



(23) 



where s — (pi +P2) 2 (pi and p 2 being the incoming four- 
momenta of the scattering particles), and 



A = 7^9 s + 9f + 4g v 



(24) 



where g s , gj and g v are respectively the scalar, fermion, 
and vector relativistic degrees of freedom. To derive this 
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FIG. 1: The graviton mode functions evaluated on the brane 
as a function of Kaluza-Klein mass m. The solid line is the 
function including the GB term, and the doted line is the one 
without GB term. The adopted parameter is x — 0.004. 



amplitude we assume that the mass of the incoming par- 
ticles is neglected. 

Going back to cosmology, the production of gravitons re- 
sults into an energy loss for ordinary matter, which can 
be expressed as 



with 



p + AHp = - dm 



C[/] = 



d 3 p„ 
(2n> 



■C[/], 



(25) 



d 3 Pi 



2 J (2n) 3 2E 1 J (2tt) 3 2E 2 



d 3 p 2 



x£|M| 2 /i/ 2 (27r) 4 <5( 4 )( Pl 



P2 - Pn 



(26) 



where /, = \ j(e Ei l T ± 1) is the Fermi/Bose distribution 
function and p m is the four-momentum of the created 
bulk graviton. The graviton production can be significant 
at high energies. So, heavy gravitons with m~T>r' 
mainly contribute to the energy loss (this is a good ap- 
proximation for a range of values for £ since the constraint 
from gravity experiment is I < 10 3 cV _1 ). Accordingly 
we have to look into the behavior of the mode function 
u m for m ^> t~ l . In the RS case we have |w m (0)| 2 = 
const for m t -1 . However, in the GB case the mode 
functions exhibit a rather nontrivial dependence on m for 
m 3> i^ 1 as shown Figure 1. For the modes m 3> £~ x , 
cqs. (20) and (21) give us 



u m (0) 



1 + 3x + ¥x 2 + x 2 m 2 e 2 



(27) 



where we have neglected the term which is smaller than 
m- 2 i~ 2 . 
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From cqs. (17), (23), (25), (26) and (27) wc find that 







k 2 A{1-j) ^ 2nT 8 
^ 1 + n 



2 9 7T 3 



dx 



>K 2 (±&x) 



7T (1 + 3 X + fx 2 + x 2 i 2 T 2 x 2 ) 



aGB(l~7)K 2 p 3/2 



for m ~ T > 



(28) 



where K 2 is the modified Bessel function of the second 
kind, and olgb,r.s is a dimensionless constant related to 
the total number of relativistic degrees of freedom. If all 
degrees of freedom of the standard model arc relativistic, 
a GB ^ 4.56 x 1CT 4 and a RS ^ 1.54 x 1CT 3 . Wc find that 
V> is proportional top 3 / 2 , not p 2 as in the RS case at 
high energy scales 14j. Note that the transition energy 
scale depends on the bulk curvature scale £ and that this 
energy scale is lower than the RS energy scale for a wide 
range of values of £ and a. 



IV. NUMERICAL ANALYSIS 

It is useful to define the dimensionless parameters, p = 
p/X, i = t/£, H = H/£, C = C/£ 2 , and a = a/£ 2 . 
The first four variables are the same as those used in 
Rcfs. [9lllo|. Using these variables, the dynamics on the 
brane are governed by the following system of differential 
equations: 



^ + AHp = 
at 



dC_ 

di 

dH _ 
di 
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4d r 4 6V 3 
_ 7(3 - \ ! 
r 4 6V 3 351/3 



4a 

2C (3 - 7)j 
36i/ 3 
I-6 1 / 3, 



(29) 
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/ - 1 - b 1 / 3 
T)(p-3) \ H 2 + 



I \ 1 + 8a [H 
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(31) 



where 



/ U 1 
V = -jV 

2a 



aj?.s(3-7)(l-7)p- 

i+3x+¥x 2 



M)p_1 forT 2 >l/x 2 ^ 2 

x (.32) 
for r 2 <T 2 «l/ X 2 ^ 2 
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FIG. 2: Evolution of e D with £ = lGeV . The 
curves from top to bottom correspond to the cases a — 



10" 



10" 9 ,10 _7 and 10" 5 . 



between tjj oc p 3 / 2 and tjj 
from above. 



We can see that the transition 
2 occurs in the second curve 



oc p 



In order to derive above equations we use eqs. (7), (14), 
(15), (16), (17) and (28). Unfortunately, it is very diffi- 
cult to find a general analytic solution to these equations, 
such as that found by Leeper et. al. for the RS case 11 1. 
So, we use the approximation of eq. (32) and solve the 
above coupled system numerically. Here, the algebraic 
constraint from the generalized Friedmann equation, 



H 2 = C + +C - 



8a 



{[& + 2(3- 7 ) 2 d(l + p) 
(3- 7 )(l + p)]} 



±V2a 



2 

2/3 



(33) 



is used to monitor numerical errors. Results from a nu- 
merical integration of this system with a variety of initial 
conditions for a and £ are shown in Figures 2-5. The ini- 
tial value of pi is taken 10 4 except that 10 4 is larger than 
the highest energy density scale Af 4 /A (2^|. In such cases 
we take pi = M 4 /A. 

Figure 2 shows the effect of increasing a while keeing 
£ fixed. Here, we define en as the ratio of dark radiation 
to standard radiation energy density: 



9C 



2(3 - 7 ) 2 r 4 p 



(34) 



The first thing to notice is that the larger a is the smaller 
e d ■ This is because the interaction between brane matter 
and the bulk gravitons weakens due to the GB term as 
can be seen in Figure 1 and the brane emits less gravitons. 
The evolution of en with d = 10 -7 is shown in Figure 
3. There is also a marked effect on e^. The increase 
of £ leads to an extension of the ip oc p 3 / 2 regime and a 
suppression of graviton emission by eq. (28). 
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FIG. 3: Evolution of e D with a = 10. The 
curves from top to bottom correspond to the cases I — 
10" 7 , 10" 5 , 1, 10" 5 , and 10 10 GeV _1 . We can also see that the 
transition between ip oc p 3//2 and ip oc p 2 occurs in the second 
curve from above. 



We have performed numerical integration of the system of 
differential equations in section IV. The different feature 
from a RS radiating brane is that the asymptotic value 

a 
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FIG. 4: The asymptotic values of cd with £ = lGeV . 



At low energies, dark radiation is produced at a negli- 
gible rate so that there is an asymptotic constant value 
for ££i as shown in Figures 2 and 3. These asymptotic 
values of ££> are shown in Figures 4 and 5. We find that 
there is an upper bound on cd and that this upper bound 
value is the final value of the RS case M : 



C GB 



< e 



RS 



3(Xrs- 



(35) 



This quantity is constrained by cosmological observa- 
tions. The number of additional relativistic degrees of 
freedom is usually measured in units of extra neutrino 
species AN V . A typical bound AiV„ < 1 [lj^ implies 
£d < 0.35. This bounds is just above the estimated 
value for the RS case. Including the GB term can help 
reduce the final value of the dark radiation term and 
hence we can easily satisfy this bound. 
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FIG. 5: The asymptotic values of eo with a — 10 7 . 



V. CONCLUSIONS 

In this paper, we have considered a GB brane that 
emits gravitons at early times, using a generalized Ads- 
Vaidya spacetime approximation. We have derived the 
dynamical equations governing the evolution of the en- 
ergy density p, the scale factor r, and the dark radiation 
parameter C in section II. In section III we have derived 
the production rate of bulk gravitons. 



for the dark radiation depends on the curvature scale I 
and the GB coupling term a. And we have demonstrated 
that the late-time dark radiation is generally suppressed 
and so cosmological limits can be easily satisfied when 
there is a GB term. 
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